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Abstract. For a given spatial distribution of the lenses and dis¬ 
tribution of the transverse velocity of the lens relative to the 
line-of-sight, a probability distribution for the lens mass for 
a single observed event is derived. In addition, similar prob¬ 
ability distributions are derived for the Einstein radius and the 
separation of the lens objects and their rotation period for a bi¬ 
nary lens. These probability distributions are distinct from the 
distributions for the lens population, as inve stigat ed e.g. by the 
mass moment method of De Riijula et al. ( 1991 ). It is shown 
that the expectation value for the mass of a certain event as de¬ 
rived in this paper coincides with the estimated average mass of 
the underlying mass spectrum as found with the mass moment 
method when only one event is considered. The special cases 
of a Maxwellian velocity distribution and of a constant velocity 
are discussed in detail. For a rudimentary model of the Galactic 
halo, the probability distributions are shown and the relations 
between the expectation values of the physical quantities and 
the event timescale are given. For this model it is shown that 
within a 95.4 %-interval around the expectation value the mass 
varies by a factor of 800. For the observed events towards the 
EMC — including the binary lens models for MACHO FMC#1 
(D ominik & Hirshfeld ^99^ and MACHO FMC#9 (Bennett et 


1996) — the results are shown explicitly. I discuss what 


al. 

information can be extracted and how additional information 
from the ongoing microlensing observations influences the re¬ 
sults. 


Key words: gravitational lensing — dark matter — Stars: low- 
mass, brown dwarfs — Galaxy: halo — planetary systems 


1. Introduction 

Some attempts have been made to obtain information about 
the mass of the lens from observed microlensing events. Un¬ 
fortunately, the mass cannot be inferred directly. Instead, the 
only relevant information directly available from a fit of a light 
curve to the observed data points is the timescale and the 
mass depends on this timescale (M cx t^) as well as on the 
position of the lens and the transverse velocity of the lens rel¬ 
ative to the line-of-sight. Since the latter parameters are both 
not observable (except for extraordinary cases where some ad¬ 


ditional information can be obtained), one can only obtain sta¬ 
tistical information on the lens mass assuming distributions of 
the lens position and the transverse lens velocity. Griest (1991), 
cited as GRI in the following, argued that the most likely mass 
distribution is that which yields the best fit to the distr ibutio n 
of timescales of the observed events. De Riljul a et al . ( 1991 ), 
cite d as R JM in the following, Jetzer & Masso ( 1994 ) and Jet- 
zer ( 1994 ) have shown that one can extract statistical moments 
of the lens mass distribution from the moments of the distribu¬ 
tion of the timescales. 

While these attempts use the distribution of timescales to 
obtain information about the mass spectrum of the lenses as re¬ 
alized in nature, I will discuss another topic in this paper. For a 
given observed event, I investigate the probability distributions 
of the lens mass and other physical quantities like the Einstein 
radius and the separation and rotation period of the lens ob¬ 
jects for a binary lens. These probability distributions give the 
answer to the question, how probable certain ranges of the con¬ 
sidered physical quantity are for the given lens system having 
produced the event. Note that this question is not answered by 
applying the methods of GRI or RJM. It is however of special 
importance for planetary systems. Also note that additional pa¬ 
rameters will enter the calculation of the event rate if one con¬ 
siders events which deviate from the point-source-point-mass- 
lens model. A binary lens will give serious problems for deter¬ 
mining the mass spectrum from the timescale distribution since 
one measures the agglomeration of two objects each from the 
mass spectrum for one event rather than a single object from 
the mass spectrum. In contrast to this, the probability distribu¬ 
tions presented in this paper also give meaningful results for 
‘anomalous’ events. 

This paper is organized as follows. In Sect. 2, it is shown 
which information directly results from a fit of a galactic mi¬ 
crolensing light curve. In Sect. 3, the relation between the 
event rate and the mass spectrum, the spatial distribution of 
the lenses, and the distribution of the relative velocity is given. 
Section 4 shows how the probability distribution for the lens 
mass can be derived. Section 5 gives results for the moments 
of further physical quantities and the probability distribution of 
these quantities around their expectation value. In Sect. 6, it is 
shown that the expectation value for the mass coincides with 
the value obtained by applying the mass moment method of 
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RJM to one event. In Sect. 7, two special forms of the velocity 
distribution are discussed in more detail: a Maxwellian distri¬ 
bution and a fixed velocity. In Sect. 8, the expectation values 
and the probability distributions for the lens mass, the Einstein 
radius, and the separation and the rotation period for binary 
lenses are shown for a simple model of the galactic halo. In 
addition, intervals corresponding to probabilites of 68.3 % and 
95.4 % are given for these quantities. In Sect. 9, the implica¬ 
tions for the observ ed eve nts towards the EMC (Alco ck et al. 

Auborg et al. |l993t D omini k & Hirshfeld [l99^ Alcock 


et al. 1996, 1997; Pratt et al. 1996 Bennett et al. 1996) are dis 


cussed explicitly. I discuss what information can be extracted 
and how the observation of more events influences the results. 


2. Information from a fit of a galactic microlensing event 


Erom a fit of the light curve to the observed data of a galactic 
microlensing event, the only dimensional parameters are the 
point of time fmax when the event occurs and the character¬ 
istic timescale related to its duration. While the point of 
time fmax does not yield any relevant information, all phys¬ 
ical quantities related to the observed event which involve a 
dimension depend on This is true not only for the ‘stan¬ 
dard model’ of Galactic microlensing — a point-mass lens 
and a point source —, but also for ‘anomalous’ events (binary 
sources and lenses, parallax effects, blending, finite size of the 
source and the lens) j] The geometry of the microlensing events 
depends on a length scale which can be chosen as the Einstein 
radius te of the lens of mass M at a distance from the ob¬ 
server, where the source is at a distance Dg from the observer 
and at a distance D^s from the lens. The Einstein radius te is 
then given by 


rE = 


AGM D^D^s 


( 1 ) 


With fi = M/Mq and x = D^/Dg, te can be written as 

te = ro \/fix{l - x) (2) 

where 


ro = 


AGMQDg 


( 3 ) 


The characteristic time scale Ie is given by 



( 4 ) 


where vx. is the transverse velocity of the lens relative to the 
line-of-sight source-observer. Note that motions of the source 
and the observer are also absorbed into this quantity. 

Eor a given Ie, the mass fi therefore follows as 


4 

A'= 72 




C" 


rQ x{l — x) Tq x{\ — x) ’ 


( 5 ) 


* Since these ‘anomalous’ events correspond to a more general 
model and the point-source-point-mass-lens model is a special case or 
an approximation, the ‘anomalous’ events are quite the normal thing! 


where is a characteristic velocity and C = v±/vc. One sees 
that /i depends on the timescale tE as well as on x and (. By 
assuming distributions of x and (( (where the distribution of (( 
may depend on x), it should in principle be possible to derive 
a probability distribution for /i. Eor doing this, I have a look at 
the event rate in the next section. 


3. Event rate and mass spectrum 


Consider a coordinate system where the lens is at rest and 
let the source move on a straight line projected onto t he len s 
plane with velocity v±. Eollowing Mao & Paczyriski ( 1991 ), 
the characteristic width w is then defined as the range of im¬ 
pact parameters for which a microlensing event occurs. Clearly, 
the width w is proportional to the Einstein radius te, so that 
w = worE{x). The event rate E is given by the product of the 
area number density of the lenses, the perpendicular velocity 
and the characteristic width of the considered type of event: 


r = nv±w. 


(6) 


Eor variable lens position, the area number density of the 
lenses has to be replaced by an integral of the volume number 
density n over the line-of-sight direction x. Eor a general lens 
population, the number density depends also on the mass fi of 
the considered objects, so that one gets 


D. 

J d^i 


( 7 ) 


as area number density of the lenses. If the mass spectrum does 
not depend on x, one can separate the x and /r-dependence by 


dn{x,n) .dnoin) 

= H{x) ■ 


d^jL 


dfj, 


(8) 


where the function H {x) follows the volume mass density p{x) 
as p{x) = poH{x), so that p{x) = po at the reference distance 
where H{x) = 1. 

The total volume number density of lenses at the reference 
distance is 


f dno{p) 

t^o = y dp , 


( 9 ) 


so that the probability for a mass p in the interval [p, p -f dp] is 

( 10 ) 


1 dnoip) 

ojip) dp = - - - dp 

no dp 


which gives the mass spectrum. 

With H{v±) dv± being the probability of finding the per¬ 
pendicular velocity in the interval [u_l , v± + dv±], one obtains 
for the event rate 

T = DgWo [ rE{x)H{x)v±H{v±) — dpdv± dx (11) 


or 


r = 


DgWo 

<M> 


dp 


r-E{x)p{x)v±H{v±)uj{p)dpdv±dx (12) 
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with the average mass <M>= 

Let Vc be a characteristic velocity and ^ . The proba¬ 

bility density for ( is then given by 

K{0=HiCv,)vc. (13) 

Note that K may depend on x. For any x, K is normalized as 

J K{C;x)dC = i. (14) 

With these definitions and rl = rg/ttat(l — x), one gets 

r = DsWoroVc J \/fix{l-x)H{x)(K{Q ■ 

(15) 

dfj, 

4. The probability density for the mass 

Eq. gives the total event rate which includes events with 
all possible timescales from the mass spectrum and the distribu¬ 
tion of the lens position and velocity. By adding an integration 
over fs and a (5-function one gets 


of the timescales for different masses /r'. By exchanging the 
roles of /r' and <e one obtains the contribution of masses in the 
interval [/i', p' + d/r'] for events with tE to the event rate as (fi' 
is called /i in the following) 

^ = DsWoroVc a J ^x{l - x) H{x) ■ 

so that i ^ gives the probability density for the mass p. 

The normalization factor is obtained by integration over p, 
which gives 

r = DsWotlP+\^c^+^ro^Pa ■ 

■ J [x{l - x)]-P H{x) KiOdxdC 
= DsWo vlP+‘^ rfPa T{-p, 2p+2), (22) 

where 

T(r, s)= J [a;(l - a;)]’' H{x) C ^(C) dC dx . (23) 


r = DsWo To Vc J \/px{l - x) H{x) C-^(C)— ■ 


For the case that the velocity distribution does not depend on 
X, the function T{r, s) separates as 


■5 



■^^/px{l- x)'] 

wC / 


dtEdfid^dx. (16) 


The event rate contribution for timescales fs in the interval 
[tE, tE + dfs] is given by dtE- 

Let us now compare different mass spectra which have only 
the mass p!, i.e. 

oc (5(p - p') • (17) 


T{r, s) = E{r)W{s), 

where 

HM =/[x(l-x)r 

lL(s) = J CKiOdC. 

The probability density for p follows as 


(24) 


(25) 

(26) 


If one assigns the same probability to any mass a-priori, one 
has uj{p) = 6{p — p'), i.e. 


dnp . , _ po 
dp ^ pMq 


S{p- p). 


(18) 


p 


More generally, one can use any explicit form of the mass 
spectrum, e.g. a power law for by using a weighting factor 

a p" 

dno 
dp 


I.e. 


{p) =apPS{p - p'), 


(19) 


1 dr 

r dp 


ro 


2p+l 


M' 


p+l/2 , 


r(-p, 2p-|-2) \tEVc 
■^Jx{l-x)H{x)CK{05(p- \ dxdC. (27) 


Note that the width wp has cancelled out. This is due to the 
fact that the fit parameters are kept fixed and only the unknown 
quantities p, x, and v± are varied. Implicitly, the same proba¬ 
bility is assigned to each parameter for the different values of 
p, X, and (. 


so that the case above corresponds to p = —1. 

For the power-law mass spectra, one obtains 

= DsWo rpVca J v^a;(l - x)H{x) ■ 

■ C K{() 5 ^E - dC,dx. (20) 

For a given p' , the probability for a timescale in the interval 
[^E, is given by dtE as a function of tE- This 

fact has been used by GRI and RJM to compare the distribution 


5. Moments of the probability distributions for physical 
quantities 

The expectation value <p> follows as 


f ^ dT ^ 

= J 

= T(-p 2p + 2) j ■ 

■C^P+^ K{C)dCdx, (28) 
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and for a general quantity G = C?(tE, a;, C) one obtains 


1 


T(-p,2p+2) 

• KiC) dC dx . 


J G(fE,a;,C) [a;(l - at)] ^ H{x) 


This means that one averages over x with the density func¬ 
tion [a;(l — x)]~'P H{x) and over C with the density function 

C 2 p+ 2 x(C). 

For the quantity G being of the form 

G(fE,a;,C) = Go(fE)[x(l-a:)]'=C' (30) 

one obtains 

<G>—Go—— — GoF{p,k,l). (31) 

T{-p,2p+2) 

For the mass, one has 


p{tE,X,C) = 


tW C 


Tg a;(l — x) ’ 


(32) 


(33) 


(34) 


so that k = —1 and I = 2 , and 

4 T(-l -p,2p + 4) 

^ T(-p,2p+2) ■ 

Forp = — 1, this gives 

_ f|tt?T( 0 , 2 ) 

^ S(l) ’ 

which is the same value as for the average mass in the mass 
spectrum (if one uses only the information from a single event), 
which can obtained by the method of mass moments described 
in Sect. ^ Note however that the probability density for p for 
a specific event and the mass spectrum are different quantities 
and that the higher moments are different. 

If the velocity distribution does not depend on x, one gets 
for the expectation value of the mass 


_ t|»? H( 0 )ir( 2 ) 

^ rS H(l) ^ 

S(0) is related to the surface mass density S by 


(35) 


Ds 


E = 


j p{D^) dDd = DsPo J H{x) dx = D^po ^( 0 ), (36) 


0 0 
and S(l) is related to the optical depth r by 

D, 


AttGD 


T = 


p{Di) dDd 


4G 2 f , ^ , 

= pqtt / H{x)x{l — x)dx 


2 DsPo , 


Using these results, and noting that 
<ul>= vl VF(2), 


(38) 


the expectation value of the mass can be written in the conve¬ 
nient form 


(29) 


<.>=*! 


(39) 


using the optical depth r, the area number density of the lenses 
E, the average square of the velocity and the timescale 

t-E- Note that <p> depends only on <u5^>, not on the form of 
the velocity distribution. 

To investigate the distribution of G around <G>, I define 


KG = 


<G> 

Using 

KtE,X,C) = P 


G _ [(r(l - x)f 
F{p, k, 1) 


' 0 


C" 




a;(l — x) ’ 

one obtains for kg the probability density 


(40) 


(41) 


1 dr 

p{HG) = ^ -1 - 

1 dKG 


1 


T{-p,2p+2) \tEV. 


ro 


2p-|-l 


• ^/x(l-x)H{x) C K{C) 
■6 {kg- kg{x,C)) dxd( 
1 


J [piKGW^^^^ 


T(-p,2p+^2) 
■H{x)CK{0 ■ 


F{p,k,l) 
F{p,kJ) 1 


T{-p,2p+2)[ 

■ [a;(l - xpG-k ^( 3 .) ^- 1+2 

'^oFip,k,l) 


•d 


[Fip,k,l)] 


[a;(l — a;)]'" 

2p + 3 


2p-l + 3 


\l\T{-p,2p+2) ^ 

■ J[x{l - x)]-P-^^+^P^T H{x) 

' KGF{p,k,l)Y^'‘ 


■K 


[x(l — x)Y 


dx . 


(42) 


Note that this distribution depends only on the form of the dis¬ 
tributions of X and C and not on any physical parameters like 
Vc or po- 

The distribution of Xg = Ig kg is given by 


(37) V'(^g) = 


In 10- [F{k,p,l) ■ 10^^]' 

|;|T(-p, 2 p+ 2 ) 
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• J[x{l - x)]-P-^^+^P^T H{x) ■ 


(43) 


All moments of the probability distribution of G can be re¬ 
duced to the function T{k,j), which separates as the product 
of S(fc) and W (j) if the velocity distribution does not depend 
on X. For the n-th moment one has 

<G'”>= y i ^ G" d/r = F{p, nk, nl). (44) 

The relative deviation of G is given by 


UG _ V<G2> _ <G>2 

<G> “ <G> 


<G 2 > 

<G >2 


- 1. (45) 


spatial distribution of the mass density and the velocity distri¬ 
bution is given (RJM). 

For the moments of the mass spectrum one gets 


nk = _:__:_ 

DgWono T(fc,2 —2fc) ® 
This gives for fc = 0 

— ^ ~ 

“ D.wonov^ ’ 

which yields 

F 1 — 

“ D,wo T(0, 2) ’ 

so that 


(51) 


(52) 


(53) 


With 


F2{p,k,l) 


F{p,2k,2l) 


[F(p,fc,Z )]2 

T{2k-p, 2l + 2p + 2) 
[T{k-p,l + 2p+2)]^ 


T{-p,2p + 2), 


one obtains 

= 'jF2{p,k,l) - 1 . 


6. The method of mass moments 


pL^ 


vf T(0,2) 4^-^ 

rfT{k,2-2k) ^ 


(46) 


For k = 1, one obtains 


fE r(0, 2) tE T(0, 2) 
^r 2 T(l, 0 )“^r 2 S(l) 


(54) 


(55) 


(47) 

since r(l, 0) = S(l), which is the same value as obtained for 
the expectation value of the lens mass for p = — 1 as derived 
above, if one considers a single observed event. 


For the moments of the fE-distribution for several events, one 
obtains following RJM 


— 1 f dr _ , 

^ r J dtE ^ 

D^woroVc f I —p;- ^ ^ 

= - ^ - J VPx{^-x)H{x)CK{C) ■ 

■tE S (tE -^ fj,x(l — x) I dtE dp d( dx 

V ^'cC J 

DsWor^''~^ f ri+lduQ^p) 

=-r-7" ’ 

■ J K(OdCdx 




:ri±i rn 
-no p 2 1 


m + 1 


1-m] ,(48) 


where 


p'^ = 


no 


p^ — f ^ 


dp 


(49) 


p^ yields for the case fc = 0, = 1, and for fc = 1 the average 

mass (in units of Mq) as 


f / ^ J 1 / dno{p) J po 

p= puj(p) dp =— — - pdp= -— . (50) 

J no J dp no Mq 


One sees that the moments of the mass spectrum can be ex¬ 
tracted from the moments of the time-scale distribution if the 


7. Special forms of the velocity distribution K{C) 
7.1. Maxwell distribution 


A distribution 

H{v±) = ^ exp {-% j (56) 

vi I ^'c J 

implies 

iT(C) =2Cexp{-C^} , (57) 

and forp(KG) and ipi^c) one has from Eqs. (p^ and (p3|) 




|.| 4 -p)JT( 2 p + 2 ) 

J[x{l - x)]-P-^^+^P^T H{x) ■ 

KG F{p,k,l)\'^^‘ 


exp < - 


dx, 


iP{Xg) = 


exp < - 


[a;(l — 

2p + 4 

21nl0- [F(fc,p,0 • 10^«] ' 
|;|S(-p)FF( 2 p + 2 ) 

• J[x{l - x)]-P-^^+^P^T H{x) 

10 >G 


[a;(l — x)]'^ 


F{p, fc, 1) 


2 / 1 ' 


dx. 


(58) 


(59) 
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For p = —1 and separation of the spatial distribution and 
velocity distribution, these equations read 


= IP(I) 


j lx(l - 


1 


l-2a: 


H{x) 


fEik + l)W{l)Y^^ |_i 

•5((KG<i>(p,7))'/'-[x(l-x)]G') • 


1 s(i) ) ■ 

•d([x(l-x)]G') . 

(67) 


exp < - 


na E{k + l)W{l)Y^‘' 


[a;(l — cc)]* 


^(1) J 


,,, , 21nl0-10^ fE{k + l)Wil)\ 

= i|E,l) i E(l) ) 


dx , (60) 

2/1 


The argument of the J-function is zero for 

X 

so that 


= i ^1 ± Y^l - 4 (kg $(p, , 


• exp • 


Z|S(1) V S(1) ) 

■ j[x{l — x)Y~^~ H[x) 

E{k + l)W{l)\ 
[x{l-xW S(l) ) 


l-2x = T\/ 1-4 (kg $(p, k))^^^ . 
This yields forp(KG) 


( 68 ) 


(69) 


2/1 ' 


P(kg) = - - —K 


-1 


dx . (61) 


Examples for these distributions are shown in Sect. where 
the galactic halo is discussed. 

7.2. Fixed velocity 

A model with a fixed velocity corresponds to 

KiO = S{C - 1). (62) 

It follows that W{n) = 1 Vn and F{p,k,l) does not depend 
on I, so that 


\k\^{-p) 
H 

H 


1-4(kg$(p, k))^^^ 
i + \ll- A{kg<^{p, k)f/'" 

i ^1 - Y^l - 4 (kg 4>(p, k)f^'" 


F{p,k,l) = ^{p,k) 

with 

With Eq. (p^, p{kg) follows ; 
P[KG) = —TTT^T- —K 


(63) 


(64) 


■(70) 


Note that p{kg) does not depend on 1. 

Since there exist (real) solutions for x only for 

[KG$(p,fc)]'/"< (71) 

Kg is restricted by kg < ^G.crit for fc > 0, or by kq > ^G.crit 
for fc < 0, where 

Since <G>= Gq 4>(p, k), the critical value of G is 
Go 


Gcrit — t«G,crit <G> — 


4k ’ 


(73) 


which is a maximum for fc > 0 and a minimum for fc < 0. 
The distribution of G can be written in terms of 


G 


J [a;(l-a;)]-^’-( 3 + 2 p)i 


KG = 


kg 


= kg4'=$(p, k), 


dx 


■H{x) ■ 

■ 6 ^(kg 4’(p, — [a:(l — x)]'^^’’^ dx . (65) 


Gcrit kg ,crit 

which yields the probability density 


p{kg) = 


1 


4d P\k\^{—p) 
1 


kg 


;-i 


1 - KG 


l/k 


(74) 


(75) 




l/k 


+ 


With the substitution 


d[a;(l - a;)]G' = - [x{l - a;)]G'-i (1 - 2cc) dx , 
one gets 


+ iJ ( 1 fl- 


and for Ag = Ig kq one gets the probability density 

1 


'(/'(Ag) = 




( 66 ) 


4d~P\k\ ^{—p) 
1 


77 I - I 1 + VI - 


1 _ lO^a/k 

+ 


(76) 


(77) 


, , [4>(p,fc)]^ 243 


-1 


+ 77 ( i ( 1 - Vl - lO^G/fc 


(78) 
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Since 


Xg = Ig RG = Ig (RG KG.crit) 


= AG-lg4’(p,A:)-2fclg2, 

(79) 

■!/(Ag) is given by 


-!/(Ag) ='0(AG + lg4>(p,A:)-f 2fclg2). 

(80) 


8. Estimates and probability distributions for physical 
quantities for a simple galactic balo model 

The quantities to be discussed here are 


- the transverse velocity v_\_, 

- the Einstein radius te (which immediately yields the pro¬ 
jected distance 2 x?’e), 

- the mass M = ^Mq, 

- the rotation period T of a binary. 


For these quantities, G'(tE, x, (), k, I, and J^(—1, k, 1) ares 
shown in Table |^. The right expression for F(—1, k, 1) is valid 
if the velocity distribution does not depend on x, p denotes the 
semimajor axis in units of Einstein radii. 

For the galactic halo, a velocity distribution of 

H{v±) = ^ exp |-^ j (81) 


can be used, where = v^. 

The mass density of halo objects is modelled as 


P{r) = Po 


+ Rgc 

q,2 _|_ j.2 


(82) 


where r measures the distance from the Galactic center, i?GC 
is the distance from the sun to the Galactic center, a is a char¬ 
acteristic core radius and po is the local density at the position 
of the sun. 

With the distance parameter x, which measures the distance 
along the line-of-sight from the observer to the EMC in units 
of the total distance Ds, one obtains 


P{x) = Po 


1 -f 


1 + 


+ Rgc “ 


(83) 


where a is the angle between direction of the Galactic center 
and the direction of the EMC measured from the observer. 
With 




Bs 

Rgc ’ 


A = l + 


^GG 


B = —2 cosa , 


H(x) can be written as 


H{x) 


A 

A + Bx^s + ' 


(84) 


(85) 


Let the halo be extended to a distance of along the line- 
of-sight. With — Dh/Roc and ^ = ^h/^s, one obtains 

€ 

^( 0 ) = / H{x)dx 


and 


A 


A + Bx^s + 


dx 


Ch 


1 


A 


J A + Bx -I- 


dx 


A 2 

Ts V4A - 52 

— arctan ■ 


arctan 


B 


2 eh + g 

V4A - ^2 


V4A - B^ 


( 86 ) 


:(l) = /x(l-.)H(.)<ix 

0 

€ 

x(l — x) 


= A 


A + Bx^s + x^^l 


■ dx 


A 

W 


Ch 


X(^s - X) 

A + Bx + 3? 


dx 


A 

§ 


73 t + 7^ + R) In 


A -f B^h + Ci 


+ B^-2A 
V4A - B2 
B 


arctan 


A 

2 gh + g 

V4A - 52 


— arctan ■ 


V4A - B^ 


(87) 


Other values of S(r) can be obtained numerically. 

Selected values for S(r) and W{s) and F{—l,k,l) are 
shown in Tables ^ ^ and p], using the values used by Paczyriski 
( 11986| )^ 

Ds = 50 kpc, Rgg = 10 kpc, Xs = 5 , Xh = 5 , 

a = 82° , a = 0 . 

From these values, one obtains for the expectation values: 


<x_L> = -vTTXc = 0.886 Xc, 


<rE> = 0.107 




210 km/s J \ 1 d 


AU, 


210 km/s J \ 1 d 




( 88 ) 

(89) 

(90) 


<p> = 2.71 • 10“ 

^ Using slightly different values for Db and Rgc, and varying the 
core radius a between 0 and 8 kpc yields estimates which differ by 
about 5 %. 
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Table 1. G{tE, x, (), k, I, and F{—1, k, 1) for v±, te, T 


G 

G{tE,xX) 

k 

i 

F{-l,k,l) 


Uc C 

0 

1 


Te 

ffi Vc C 

0 

1 


M 

T 

t|“o 

-1 

1 

2 

T(0,2) S(0)iy(2) 

7-2 X(l-X) 

^ v/p® 1 e Ds Vc x{l - a:) Vc 

1 

H(l) H(l) 

T(§.1) H(§)W(1) 

2 

2 

2(1) 2(1) 


Table 2. Selected values for ^(r) 


r H(r) 

0 0.305 

i 0.105 

1 0.0407 

1 0.0168 

2 0.00721 


Table 3. Selected values for W{s) 


s W{s) 

-i 1.225 

0 1.000 

i 0.906 

1 0.886 

1 0.919 

2 1.000 


Table 4. J^(—1, k, 1) for v±, te, ft, T 


G 

Fi-l,k,l) 

v_L,r-E 

0.886 

p 

7.49 

T 

0.374 


<T> = 


^210 km/s 

J,2 

With po = one obtains 


1/2 ,^, 1/2 

Id 


a. 


Po = 8.16-10 


-3 


210 km/s J pc^ ’ 


MfT 


© 


(91) 


(92) 


so that for Vc = 210 km/s, one gets S = 124 Mq/pc^ and 
r = 5 • 10-^. 

The distribution of = p/ <p> is given by 


p{k^) = 


( 0 ) 


(S(l))^ 


a;^(l — x)'^ H{x) 


■ exp n^x{\ — a:) ^ dx . 


(93) 


and the probability density 

t(l))2 


lO'^'^lnlO / x^{1 — x)'^ H(x) 


exp ^ 10^'‘a;(l — x) ^ dx . 


(94) 


From Eq. ( |93[ ) the probability density P{p) for the mass p (in 
units of Mq) follows as 

^ ^ ^ ~ ' 
p Tq a;(l — x) 


■ exp 




dx , 


(95) 


which differs from the probability given by Jetzer & Masso 
(1994, Eq. (8)), by a factor /i^/f|. Note that this probability 


(96) 


density has to be of the form 

''E \ ''E 

since ft cx ^e, to ensure normalization for any Ie- 

Eor the other quantities estimated the probability densi¬ 
ties are given by Eqs. (^) and (^). The probability densities 
p(kg) and '0 (Ag) for te, ft, T are shown in Eigs. |i| to |[ Note 
that Te follows the velocity distribution, because fc = 0. In 
the diagrams for ij}{Xc), symmetric intervals around <G> are 
shown which give a probability of 68.3 % and 95.4 % respec¬ 
tively. The bounds of these intervals are also shown in Table 
The bounds are much larger for te and again much larger 
for ft than for T, which is due to the wide distribution of the 
velocity and p cx ^e cx C, while T cx v^. The smallest 
and the largest value in the 95.4 %-interval differ by a factor of 
about 800 for ft, 16 for te and 5 for T. 


9. Application to observed events 

In this section I show the application of the method described 
here to the observed events towards the EMC. The first events 
have been claimed by EROS (Aubourg et al. 

EROS#l and #2, and MACHO (Alcock et al. 


19931, namely 


1993), namely 


MACHO LMC#1, in 1993. The fit with a point-mass lens and 
point source for MACHO LMC#1 showed a discrepancy near 
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Table 5. The bounds of symmetric intervals on a logarithmic scale around <G> which correspond to probabilities of 68.3 % and 
95.4 % 


G 

AA68.3 

AA95,4 

■^Q — ^^68.3 

10^^68.3 

10~'^^95.4 

^qAA95.4 

v±, Te 

0.2428 

0.6111 

0.572 

1.75 

0.244 

4.09 

p 

0.5900 

1.454 

0.257 

3.89 

0.0351 

28.5 

T 

0.1588 

0.3719 

0.694 

1.44 

0.425 

2.35 


Fig. 1. The probability density ) = p{K,r-E,) for= v±/<v±>= te/ <rE>= Krs (left) and The probability density 
) = V’('^rE) with symmetric 68.3 % and 95.4 % intervals around <r;_L> or <rE> (right) 


Fig.2. The probability density for = p,/ <p> (left) and the probability density with symmetric 68.3 % and 

95.4 % intervals around <p> (right) 


the peak which has been solve d with mode ls involving a binary 
lens by Dominik & Hirshfeld ( |l994l |l996| ). The MACHO col¬ 
laboration had found t wo ot her candidates, MACHO LMC#2 
and #3, (Alcock et al. 1996| ) which have been meanwhile dis¬ 
missed. In addition, they have claimed the existence o f 7 ad - 
ditional events, MACHO LMC#4 to #10, (Pratt et al. |1996|), 


where MACHO LMC#9 is due to a binary lens (Bennett et al. 


1996). In the MACHO data taken from 1996 to March 1997, 5 


additional LMC candidates showed up (Stubbs et al. |1997[ ). 

It has been shown that the E ROS# 2 event involves a pe¬ 
riodic variable star (Ansari et al. 1995) and that EROS#l in¬ 


volves an emission line Be type star (Beaulieu et al. 1995), so 
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Fig. 3. The probability density p{kt) for kt = T / <T> (left) and the probability density '(/'(At) with symmetric 68.3 % and 
95.4 % intervals around <T> (right) 


that both EROS#l and EROS#2 involve a rare type of stars 
which makes t hese e vents suspect as microlensing candidates 
(e.g. Paczyhski |l996| ). In addition, the MA CHO LMC#10 event 


is likely to be a binary star (Pratt et al. 
199 ^. 


; Alcock et al. 


By assuming that the lens is in the galactic halo and using 
the halo model of the last section, expectation values for the 
desired quantities can be obtained by inserting the fit parame¬ 
ters into Eqs. ( ^8| ) to (^l]). Table ^ shows the expectation val¬ 
ues for the Einstein radius and the mass for the events EROS#l 
and EROS#2, MACHO LMC#4... #8 and #10, whereas the re¬ 
sults for the binary lens events MACHO LMC#1 and MACHO 
LMC#9 are shown in Table 0. Eor MACHO LMC#1 six differ- 
ent binary lens models are shown (Dominik & Hirshfeld 1996; 
Dominik [l9^ . Note that the lens for MACH O LMC#9 prob¬ 
ably resides in the EMC (Bennett et al. 19961. 


The parametrization used is that of Dominik & Hirshfeld 


(19961: Thd denotes the projected half-distance between the 


lens objects in the lens plane, where rud = 2xr’E, rni denotes 


the mass fraction in lens object 1. 


.( 2 ) 


radius corresponding to the mass of object 2 and ’ the char- 


denotes the Einstein 

( 2 ) 


acteristic time corresponding to r 


.( 2 ) 

'E 


( 2 ). 


' = r-E Vl - mi , 


l( 2 ) 


(97) 

(98) 


i-E = t-E VI - mi. 

Similarly x denotes the projected half-separation in units of 
Einstein radii corresponding to the mass of object 1, so that 


^(i) = 




X- 


(99) 


pi and p ,2 denote the mass in units of Mq for objects 1 and 2 
and the mass ratio r is given by 


1 — mi 

mi 


(100) 


Table 8. MACHO LMC#1: Estimates of physical parameters 
of the binary lens models using a mass ratio r at the upper 2-cr- 
bound of r and corresponding values of and % . 



BA 

BAl 

BA2 

BA3 

is/d 

57.67 

86.93 

29.94 

228.6 

4^Vd 

17.63 

15.20 

17.61 

16.49 

r 

0.1031 

0.0315 

0.529 

0.00523 

X 

2.291 

2.229 

1.735 

2.222 

x« 

2.406 

2.264 

2.146 

2.228 

<rE> /AU 

6.20 

9.34 

3.22 

24.6 

<4">> /AU 

1.90 

1.63 

1.89 

1.77 

<2rhd> /AU 

28.4 

41.7 

11.2 

109 

<p> 

0.90 

2.0 

0.24 

14 

<pi> 

0.82 

2.0 

0.16 

14 

<P2> 

0.084 

0.063 

0.084 

0.074 

<Tmin> /a 

69 

82 

33 

132 


Since the true semimajor axis a = p te is not yielded by the 
fit, T is estimated using p = which corresponds to a mini¬ 
mal value Tmin, because p > x for gravitationally bound 
system and T cx 

The distribution of the physical quantities as well as sym¬ 
metric intervals around the expectation value with probabilities 
of 68.3 % and 95.4 % are shown in the previous section. 

The model BA3 has previously been omitted (Dominik & 
Hirshfeld 199^ due to the fact that the mass ratio between the 
lens objects is very extreme (4 • 10“^^), which corresponds to 
unphysical values (see Table However, the uncertainty of 
the mass ratio is very large for extreme mass ratios (Dominik 
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Table 6. Expectation values for the point-source-point-mass-lens events towards the EMC. M # denotes MACHO events, while 
E # denotes EROS events. 



M#4 

M#5 

M#6 

M#7 

M#8 

M#10 

E#1 

E#2 

tfi/d 

23 

41 

44 

58 

31 

21 

27 

30 

<rE> /AU 

2.5 

4.4 

4.7 

6.2 

3.3 

2.2 

2.9 

3.2 

<At> 

0.14 

0.46 

0.52 

0.91 

0.26 

0.12 

0.20 

0.24 


Table 7. Expectation values of the physical parameters and used fit parameters for the 6 binary lens models for MACHO LMC#1, 
denoted by BE, BLl, BA, BAl, BA2, and BAS, and for MACHO LMC#9. 



BE 

BLl 

BA 

BAl 

BA2 

BA3 

LMC#9 

tfi/d 

16.27 

17.53 

685 

155 

35.7 

2.62 • 10^^ 

143.4 

4^Vd 

— 

— 

17.57 

15.15 

17.72 

16.36 

— 

mi 

0.463 

0.557 

0.99934 

0.9904 

0.75 

1 - 4 • 10“^® 

0.620 

r 

1.16 

0.795 

6.6 ■ 10"'^ 

9.7- 10"® 

0.33 

3.9 • 10“^® 

0.613 

X 

0.20 

0.22 

2.41 

2.21 

1.82 

2.24 

0.83 


— 

— 

2.41 

2.22 

2.10 

2.24 

— 

<rE> /AU 

1.75 

1.88 

73.6 

16.7 

3.84 

2.82 • 10“ 

15.3 

/AU 

— 

— 

1.89 

1.63 

1.90 

1.76 

— 

<2rhd> /AU 

0.71 

0.82 

355 

73.6 

14.0 

1.2 • 10^^ 

2.54 

<At> 

0.072 

0.083 

127 

6.5 

0.34 

1.9 • 10^^ 

5.6 

<Ati> 

0.033 

0.046 

127 

6.4 

0.26 

1.9 • 10^^ 

3.5 

<A‘2> 

0.038 

0.037 

0.084 

0.062 

0.086 

0.073 

2.1 

<Tmin> /a 

0.98 

1.1 

257 

107 

39 

1.43 • 10’’ 

24 


& Hirshfeld 1996 ) due to the fact that the lens behaves nearly 
like a Chang-Refsdal lens. The same degenaracy has recently 
been rediscovered in the contex t of lensing by a star with a 
planet by Gaudi & Gould (1997). Table ^ shows the estimates 
for the 4 wide binary lens models (BA, BAl, BA2, BAS) where 
ht parameters at the upper 2-cr-bounds of the mass ratio r have 
been used. It can be seen that the expectation values for the 
separation and the mass change dramatically for the fits with 
small values of r, leaving much room for speculations about 
the nature of the lens object. 


Recently, Rhie & Bennett (1996) have speculated about a 
planetary companion in the MACHO LMC#1 event. However, 
as shown in this and the last section, there are two fundamen¬ 
tal uncertainties (beyond the fact that there are only a few data 
points near the peak, which could have been solved by a denser 
sampling), namely the uncertainty in the mass ratio and the un¬ 
certainty due to the unknown lens position and its velocity. The 


planetary model of Rhie & Bennett (1996) corresponds to my 


model BAl, where the expectation value of the mass of the low- 
mass object is about 0.06 Mq. Taking into account an uncer¬ 
tainty of a factor of 30, a value of 2 Jupiter masses is reached 
just at the 2-cr-bound. As shown in Table |^, the expectation 
value for the mass of the low-mass object is about the same 


for all of my wide binary lens fits, so that there are the same 
prospects for a planet as the low-mass object for all of these 
models. However, the high-mass object will be different. 

The estimates of physical quantities like the mass of the 
lens object(s), their separation and the rotation periods along 
with the uncertainties involved are needed to reveal the physi¬ 
cal nature of the lens for each observed event. One can check 
whether the mass range is consistent with the assumption of a 
dark object, and of which nature the lens should be in this case 
(brown dwarf, white dwarf, neutron star, ...). The determina¬ 
tion of the mass range is crucial for claiming the existence of a 
planetary companion. 

One can also check whether it is consistent to use a static 
binary lens m odel rather than one with rotating binary lens (Do¬ 
minik 1997 ). Moreover, by comparing estimates for different 
lens populations (e.g. the galactic halo and the EMC halo) one 
may obtain indications to which population the lens belongs. 

In contrast to estimates on the mass spectrum, there is no 
direct influence from the ensemble of observed events on the 
estimates for a specific event. However, there is an indirect in¬ 
fluence, since the ensemble of events gives information about 
the mass spectrum, which in turn can be used to get a more 
accurate estimate for each observed event. However, a lot of 
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events (~ 100) are needed to get accurate information on the 
mass spectrum (Mao & Paczyhski |1996| ). To be used for the es¬ 
timates for a specific event, some of the higher mass moments 
should have been determined. Note that there will remain large 
uncertainties in the mass of a specific event even if the mass 
spectrum is known (unless it contains a sharp peak), since due 
to the broad distribution of the velocity, the range of timescales 
t'E for a certain mass is broad, so that the mass range for a spe¬ 
cific event in turn may also be large. 

While a white dwarf sce nario is preferred by the LMC ob¬ 
servations (e.g. Prat t et a l. 1996 ), a brown-dwarf scenario is 
not ruled out (Spiro 1997 ) if one considers a flat halo. Though 
there are some restrictions on the average lens mass from the 
observed events, which will improve with more data, it will re¬ 
main highly uncertain into which mass regime a certain event 
falls. 

Since the ongoing microlensing observations constrain the 
galactic models and therefore give rise to more accurate deter¬ 
minations of the structure of the lens populations, the knowl¬ 
edge on specific events will also be improved by this. 


Acknowledgements. 1 would like to thank the MACHO collaboration 
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Hirshfeld for a critical reading of the manuscript, and P. Jetzer for 
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